Entanglement-screening by nonlinear resonances 
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We show that nonlinear resonances in a classically mixed phase space allow to define generic, 
strongly entangled multi-partite quantum states. The robustness of their multipartite entanglement 
increases with the particle number, i.e. in the semiclassical limit, for those classes of diffusive noise 
which assist the quantum-classical transition. 
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The ultimate success of the quantum information and 
computation program will depend on our theoretical un- 
derstanding and experimental control of quantum entan- 
glement. In order to compete with the available classical 
supercomputing resources, a future quantum computer 
will have to be composed of a large, i.e., at least meso- 
scopic number of qubits, and their coherence will need to 
be preserved over a significant period of time. Contem- 
plating the fact that entanglement is a manifestation of 
multi-particle coherence, and that the density of states 
explodes exponentially with the particle number, it is 
easy to appreciate the dimension of the challenge ahead. 

So far, little is known about entanglement in open 
quantum systems - where "open" refers to the unavoid- 
able coupling to uncontrolled degrees of freedom in the 
"environment" , which eventually induces decoherence in 
the system dynamics. Only recently [l|, [3, H, 0, H, S, 13 
have there been first theoretical and experimental at- 
tempts to characterize entanglement dynamics under de- 
coherence, but a sufficiently general picture still has to 
emerge. In particular, most studies did so far focus on 
specific classes of highly entangled W, GHZ, or cluster 
states, and on their specific robustness against certain 
sources of decoherence [1, . As the number of particles 
increases, the faithful experimental generation and prob- 
ing of these states tends to become more and more diffi- 
cult, with rapidly increasing experimental overhead. Fur- 
thermore, whether and in which sense their entanglement 
properties can be considered as "generic" is a largely open 
issue, given the complicated topology of state space. 

In our present contribution, we will adopt a different 
perspective, which imports some generic features of quan- 
tum dynamics with underlying mixed, regular-chaotic 
phase space structure, i.e., from quantum chaos In 
contrast to earlier studies of the impact of mixed phase 
space dynamics I 111 . Il2l [13| and nonlinear light-matter 
on bipartite entanglement, we are 
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interaction 

here interested in the multipartite limit of large particle 
numbers. This is of crucial importance in the context of 
entanglement scaling alluded to above, and will be iden- 
tified with the semiclassical limit of progressively finer 



(quantum) resolution of classical phase space structures, 
by a suitable definition of many-particle basis states in 
terms of classical phase space coordinates. We will see 
that nonlinear resonances, which are ubiquitous in clas- 
sical Hamiltonian systems [l6j . naturally define strongly 
- though not maximally ~ entangled multipartite quan- 
tum states. This non-optimality is compensated by the 
nonlinear resonance structure providing a natural shel- 
ter against certain types of decoherence: Indeed, the ro- 
bustness of the associated multipartite entangled states 
is found to increase with the number of particles, i.e., in 
the semiclassical limit. 

We start out with a system of k qubits that lives on 
a Hilbert space with tensor structure H. — Tii ® 0.2 ® 
. . . ® Tik, where each factor space TLj has dimension 
two. The "computational basis" {|«)}, which spans 7Y, 
is given by the /c-particle product states (i.e., by binary 
fc-strings, which we identify with the binary represen- 
tations oi i = Q...N), and has dimension N = 2^ . 
Maps on classical phase space can be implemented ef- 
ficiently in such systems after suitable identification of 
the computational with the position basis {l^i)}, upon 
the substitution = i 0, H IH • Then, the corre- 
sponding momentum basis {JPi)} is given by the discrete 
Fourier transform of [20[. Thus, the phase space is 

wrapped on a torus, with a phase space area 1/A'^ occu- 
pied by a single position basis state. Coherent (or mini- 
mum uncertainty) states are obtained as Gaussian wave 
packets of width both in position and momentum 

j2l| . Accordingly, the effective Planck constant is given 
by hgff = 1/2ttN, and the semiclassical limit /ieff 
of arbitrary phase space resolution by a single quantum 
state is approached as the particle number k tends to 
infinity. 

Given this quantum coarse graining of phase space, we 
will now monitor the time evolution of the multipartite 
entanglement of a quantum state initially prepared as a 
minimum uncertainty Gaussian wave packet, launched at 
different positions on the torus (which we will unfold as 
a unit square, for the ease of illustration). We propagate 
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the wave packet by the unitary operator 

Jjt _ ^^iNx2Cos(2irq)^iNxicos{2-!Tp)Y j^j^^ 

(with t integer) which foUows from the quantization of 
the classical Harper map 



P' = P-Xi sin(27rq) 



9' = <1 + X2 sin(27rp' 



(modi) , 



(2) 



Equation ^ describes a particle subject to a periodic 
impulse, with a position-dependent amplitude [22l |. and 
exhibits a transition from regular to chaotic dynamics 
as the kicking strengths XI1X2 £^rc tuned from XI1X2 ^ 
0.11 to XijX2 ^ 0.63. Here, we employ the fixed values 
Xi = X2 = X = 0.4964, what defines a mixed phase space 
structure as depicted in the inset of Fig. [TJ 




FIG. 1: (Color online) Evolution of multipartite concurrence 
Ck for two different numbers of qubits: k = 5 (top), k = 8 
(bottom). Open symbols refer to unitary dynamics, while 
filled symbols represent the evolution under the diffusive noise 
described by Eq. (|3]). Squares correspond to an initial condi- 
tion inside the nonlinear resonance island, triangles to initial 
conditions within the chaotic sea, in the classical phase space 
spanned by p and q (see inset). For different qubit numbers, 
the ratio of noise strength to the effective size of Planck's 
quantum, e/h^s, is kept constant (e = 0.04 for k — 5, and 
e = 0.005 for fc = 8). Black arrows indicate the value of Ck 
for fc-partite GHZ states. 

The quantum evolution ([1]) modifies the decomposi- 
tion of the evolved state \tpt) in the computational ba- 
sis, and thus the entanglement of the k degrees of free- 
dom defined by the constituent qubits. In order to assess 
the nonclassical correlations inscribed into the evolved 



state, we use the fc-partite concurrence Ck as defined in 
[l[. On pure states, this quantity is given by the square 
root of a balanced average over the squared concurrence 
of all nontrivial bipartitions of the fc-set under scrutiny, 
and vanishes exclusively for fc-separable states [H [23| . 
Moreover, it has the particularly advantageous property 

CM ® \^p)k-i) = Cfc_i(|V)fc-i), V|</>) e nu\^)k-i e 

7^2 (81 ... (81 Hk ■ The latter allows to compare the entan- 
glement inscribed into quantum states composed of an 
increasing number of subsystems. Furthermore, Ck has 
a generalization for mixed states (through the convex roof 
construction [l[), which we will make use of below. 

We start out with a short inspection of the entangle- 
ment dynamics under purely coherent dynamics, to set 
the scene. The evolution is generated by application of 
[/*, where the integer t counts the number of applications 
and defines a discrete time. The open symbols in Fig. [1] 
represent Ck{t), for two different initial positions (also in- 
dicated by filled triangles and squares in the inset) of the 
initial minimum uncertainty state in phase space - either 
within an elliptic island {q ^ p = 0.25; open squares), or 
within the chaotic phase space component [q = 0.25, 
p = 0.0; open triangles). Since the computational basis 
is encoded in position states, both initial conditions de- 
fine the same initial value Cfc(O) (6*5(0) = 1.074386 and 
6*8(0) = 1.316826, respectively) of concurrence. Since, 
on the quantum level, chaotic dynamics is tantamount of 
strong coupling in any basis, it immediately follows that 
Ck will increase rapidly for the initial condition placed 
in the chaotic domain, and saturate once equilibrated 
over the chaotic eigenstates U, 12, of the quantized 
Harper map - this is indeed observed in the figure. Note, 
however, that the saturation level does not coincide with 
the maximal possible value of Ck on the pure states, but 
rather with its most probable value (which approaches 
the maximal value in the limit k ^ ao [l^l). Also note 
that, for fc > 3, this value is larger than the fc-partite 
entanglement of GHZ states, which are maximally en- 
tangled only in the special case fc = 3 [361 . Indeed, 
Ck=5 = 1-369286, and Cg|^ = 1.408798 (indicated by 
black arrows in Fig. [1]), approaching our minimal un- 
certainty state's initial multipartite entanglement with 
increasing fc. Finally, the qualitative behaviour of Ck{t) 
for an initial state lying in the chaotic sea depends only 
weakly on the number of qubits, as evident from a com- 
parison of both panels in Fig. [TJ For fc = 8 (correspond- 
ing to a Hilbert space dimension N ~ 256), the time 
evolution is very smooth, while for fc = 5 (Hilbert space 
dimension N = 32) still some fiuctuations - essentially a 
finite size effect - are observed. 

For the state initially placed within the elliptic is- 
land, size does matter: for fc = 5, the initial coherent 
state cannot be well accomodated within the elliptic is- 
land in phase space (due to the finite size of hes), and 
exhibits non-negligible tunneling coupling to its chaotic 
environment. Consequently, as time proceeds, the co- 
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herently evolved state spreads more and more over the 
chaotic phase space component, and its entanglement fi- 
nally reaches essentially the same value as for the initial 
condition within the chaotic domain, just after consid- 
erably longer time - essentially determined by the rele- 
vant tunneling matrix elements (which, in general, will 
be strongly fluctuating under small parameter changes 



2^,1261). In contrast, for fc = 8, tunneling from the island 
into the chaotic sea occurs on a much longer time scale 
(which, on average, increases exponentially with Hes), 
and remains invisible on the time scale covered in Fig. [TJ 
The small oscillations of Ck=8{t) are due to the spreading 
of the initial wave packet along the regular island's tori. 

This screening of the initial state from the chaotic sea 
when initially placed within the elliptic island, more and 
more efficient with increasing particle number, has an 
immediate consequence for the robustness of the state's 
multipartite entanglement under the influence of deco- 
hercnce, as illustrated by the filled symbols in Fig.[TJ For 
chaotic initial conditions, an initial rise of Ck is rather 
quickly overruled by the loss of multiparticle coherence 
and hence of entanglement, and this is once again largely 
independent of fc. However, for initial conditions within 
the island, k ~ 5 again leads to asymptotically the same 
behaviour as for the chaotic initial condition, while fc = 8 
induces entanglement dynamics almost completely unaf- 
fected by the noise. Thus, for sufficiently large fc, equiv- 
alent to sufficiently small hcff, and correspondingly sup- 
pressed tunneling rates, the classical nonlinear resonance 
creates strongly entangled multipartite states which, in 
addition, are robust against noise. This is further illus- 
trated in Fig. [51 where Ck=8{t = 16) is plotted for differ- 
ent initial conditions, in the absence and in the presence 
of noise: Clearly, entanglement is robust when shielded 
by the resonance island. While chaotic dynamics produce 
slighty stronger entanglement, this is significantly more 
fragile under decoherence. 

Given the above, some remarks on the applied noise 
are in order, since any statement on the robustness of 
some kind of coherence must depend on the nature of 
the noise. Here, we used Gaussian noise as described in 
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28[ , which can be implemented experimentally by an 
enlargement of the quantum register by a set of suitably 
initialized ancilla qubits [2^ . It can be written as a map 



(3) 



qp 



which is applied once after each application of U^. 
c^(j),q) is the discrete Fourier transform of Ct{fi,h') = 



exp 



f.N^ 



(sin^ [nfi/N] 



[nv/N]) 



and is very 



close to a periodic Gaussian of width e/(27r), centered 
around {q,p) = (0,0). The Tqp are unitary translation 
operators on the torus (s^l- 

The action of Se(/5) is easily understood in phase space: 
with high probability, the state is left untouched, while 




FIG. 2; (Color online) Eight-partite concurrence Ck=8 as a 
function of initial momentum {p = 0.0, 0.05, . . . , 0.95, 1.0), af- 
ter 16 iterations of the Harper map, with x = 0.4964. Open 
symbols correspond to unitary evolution, while filled symbols 
refer to unitary evolution amended by diflusive noise, Eq. ([3}. 
Squares represent the initial position q — 0.5, and triangles 
q — 0.25. Peaks and dips of Ck=s{t) are located exactly at 
the center of the nonlinear resonance island in classical phase 
space. 



with weight Ce(p, q), locally in phase space, every possible 
translation is generated. This noise is similar to a high 
temperature bath of oscillators producing both diffusion 
and decoherence, and steeres the quantum dynamics into 
the semiclassical limit [3l|. In Wigner phase space rep- 
resentation, diffusion induces broadening and blurring of 
the contour of the state, while decoherence wipes out the 
interference fringes, and eventually transforms the state 
into a mixed (classical) state. Given the fc-dependence 
of heS, the strength e of the noise was scaled such as to 
keep e/fieff constant, in the above plots. The observed 
robustness of entanglement stems from the local action 
of Se(p), which respects the classical phase space struc- 
ture and leaves the initial state within the elliptic island 
effectively as a fixed point of the evolution, in the limit of 
large fc (small Ses). The dramatically different dynam- 
ics of Ck=s within and outside the island is a (multipar- 
tite) manifestation of enhanced decoherence in classically 
chaotic as opposed to regular systems (32j . for this spe- 
cific type of noise, and suggests the emergence of a sub- 
space which is shielded against disentanglement, in the 
semiclassical limit. 

In contrast, we may choose non-local noise sources such 
as defined by the multipartite phase damping channel 
(PDC), 

Sf^^(p)^(l-6)p + 65^p,,|0(z|, (4) 

i 

or the generalized depolarizing channel (DPC) 

Sr^(p) = (l-6)p + 6/. (5) 

For a computational basis identified with position states, 
both these maps can be expressed in terms of transla- 
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tions on the torus [33[, though now with equal (rather 
than Gaussian, see ^) weight. Thus, they completely 
obliterate the classical phase space structure, and cou- 
ple different eigenstates of the quantized Harper map, 
irrespectively of their localization properties in regular 
or chaotic phase space domains. Accordingly, given the 
same initial conditions as in Figs. [1] and [2l concurrence 
decreases monotonously with time, independently of the 
initial condition, except for the typical short time tran- 
sient (see also Fig.[T]) observed for the chaotic initial con- 
dition - sec Fig. m 
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FIG. 3: (Color online) Evolution of multipartite concurrence 
Cfc=8 in the presence of phase damping (triangles, (Eq. Q) 
and depolarizing noise (squares, Eq. ([5])), with e — 0.04 in 
Eqs. (|4)| and ((5}. Filled and open symbols refer to initial 
conditions in the chaotic domain and within the resonance 
island, respectively, precisely as in Fig. [T] 

To conclude, we have shown that minimum uncertainty 
states induce multipartite entanglement in the associ- 
ated computational basis, robust against the action of 
diffusive Gaussian noise, when launched within a nonlin- 
ear resonance island. They can be produced efhciently 
as ground states of the Harper Hamiltonian [l^, sub- 
sequently translated using torus translation operators. 
The latter generate modular additions with controlled 
phase shifts on the register qubits, and can be imple- 
mented, e.g. in ion trap experiments (b^. [35l|. Further- 
more, given the robust entanglement evolution for initial 
conditions within a regular island, also those eigenstates 
of the Harper Hamiltonian which are anchored to the 
classical regular island exhibit the same robustness prop- 
erties. Thus, robust multipartite entangled states can be 
defined through the resonance condition which defines 
the regular island, a ubiquitous feature of Hamiltonian 
systems with mixed classical phase space structure. 

F G.-M. was partially supported by CONICET (Ar- 
gentina) and by the EC IST-FET project EuroSQIP. 
F. M. acknowledges financial support of Alexander 
von Humboldt foundation. 



* Electronic address: lgarcia@irsamc. ups-tlse.fr] 
[1] F. Mintert, A. R. R. Carvalho, M. Kus, and A. Buchleit- 

ner, Phys. Rep. 415, 207 (2005). 
[2] K. Zyczkowski et al., Phys. Rev. A 65, 012101 (2001). 
[3] B. V. Fine et al., Phys. Rev. B 71, 153105 (2005). 
[4] A. R. R. Carvalho et al., (20 05), URL 

|http : //arxiv . org/abs/quaiit-ph/0510006] 
[5] T. Yu and J. H. Eberly, Phys. Rev. Lett. 97, 140403 
(2006). 

[6] C. Petitjean and Ph. Jacquod, Phys. Rev. Lett. 97, 
194103 (2006). 

[7] G. Puentes, D. Voigt, A. Aiello, 

and J. P. Woerdman (2006 ), URL 
|http : //arxiv . org/abs/quant-ph/0607014| 

[8] A. R. R. Carvalho, F. Mintert, and A. Buchleitner, Phys. 

Rev. Lett 93, 230501 (2004). 
[9] W. Diir and H.-J. Briegel, Phys. Rev. Lett. 92, 180403 

(2004). 

[10] J. Madronero et al, in Advances in Atomic, Molecular, 

and Optical Physics, edited by G. Rempe and M. O. 

Scully (Elsevier, Amsterdam, 2006). 
[11] K. Furuya, M. C. Nemes, and G. Q. Pellegrino, Phys. 

Rev. Lett 80, 5524 (1988). 
[12] P. A. Miller and S. Sarkar, Phys. Rev. E 60, 1542 (1999). 
[13] H. Fujisaki et al., Phys. Rev. E 67, 066201 (2003). 
[14] A. Kowalewska-Kudlaszyk and W. Leohski, Phys. Rev. 

A 73, 042318 (2006). 
[15] M. K. Olsen, Phys. Rev. A 73, 053806 (2006). 
[16] A. J. Lichtenberg and M. A. Lieberman, Regular and 

Chaotic Dynamics (Springer- Verlag, New York, 1991). 
[17] A. J. Scott and C. M. Caves, J. Phys. A: Math. Gen. 36, 

9553 (2003). 

[18] C. Miquel, J. P. Paz, and M. Saraceno, Phys. Rev. A 65, 
062309 (2002). 

[19] J. P. Paz, A. J. Roncagiia, and M. Saraceno, Phys. Rev. 

A 69, 032312 (2004). 
[20] J. H. Hannay and M. V. Berry, Physica ID, 267 (1980). 
[21] S. Nonnemacher, Nonlinearity 10, 1569 (1997). 
[22] P. Leboeuf, J. Kurchan, M. Feingold, and D. P. Arovas, 

Phys. Rev. Lett. 65, 3076 (1990). 
[23] R. Demkowicz-Dobrzanski et al., Phys. Rev. A 74, 

138611 (2006). 
[24] C. Viviescas, et. al (in preparation). 
[25] S. Tomsovic and D. UUmo, Phys. Rev. E 50, 145 (1994). 
[26] C. Ehschka and P. Schlagheck, Phys. Rev. Lett 94, 

014101 (2005). 
[27] S. Nonnenmacher, Nonhnearity 16, 1685 (2003). 
[28] L Garci'a-Mata and M. Saraceno, Phys. Rev. E 69, 

056211 (2004). 
[29] L Garci'a-Mata (Unpublished). 

[30] J. Schwinger, Proc. Natl. Acad. Sci. 46, 570 (1960). 
[31] A. R. R. Carvalho et al, Phys. Rev. E 70, 026211 (2004). 
[32] W. H. Zurek and J. P. Paz, Phys. Rev. Lett 72, 2508 
(1994). 

[33] M. L. Aolita et al., Phys. Rev. A 70, 062301 (2004). 
[34] V. Vedral et al., Phys. Rev. A 54, 147 (1996). 
[35] C. Miquel et al., Phys. Rev. A 54, 2605 (1996). 
[36] This corrects an earlier statement in [^. 



